The Peregrine rogue waves induced by interaction between the continuous wave and 

soliton 
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Based on the soliton solution on a continuous wave background for an integrable Hirota equation, 
the reduction mechanism and the characteristics of the Peregrine rogue wave in the propagation of 
femtosecond pulses of optical fiber are discussed. The results show that there exist two processes 
of the formation of the Peregrine rogue wave: one is the localized process of the continuous wave 
background, and the other is the reduction process of the periodization of the bright soliton. The 
characteristics of the Peregrine rogue wave are exhibited by strong temporal and spatial localization. 
Also, various initial excitations of the Peregrine rogue wave are performed and the results show that 
the Peregrine rogue wave can be excited by a small localized (single peak) perturbation pulse of 
the continuous wave background, even for the nonintegrable case. The numerical simulations show 
that the Peregrine rogue wave is unstable. Finally, through a realistic example, the influence of the 
self-frequency shift to the dynamics of the Peregrine rogue wave is discussed. The results show that 
in the absence of the self-frequency shift, the Peregrine rogue wave can split into several subpuslses; 
however, when the self-frequency shift is considered, the Peregrine rogue wave no longer splits and 
exhibits mainly a peak changing and an increasing evolution property of the field amplitude. 

PACS numbers: 42.65.Tg, 42.81. Dp, 42.79.Sz 



I. INTRODUCTION 

A rogue wave is an oceanic phenomenon with ampli- 
tude much higher than the average wave crests around 
them So far, this phenomenon has not been under- 
stood completely due to the difficult and restricted ob- 
servational conditions. Therefore a great deal of atten- 
tion has been paid to better understanding their physical 
mechanisms. It has been suggested that the rogue waves 
appearing in the ocean are mainly caused by wave-wave 
nonlinear interaction, such as a modulation instability 
of the Benjamin- Feir type @-@|. Recently, rogue waves 
have been observed in optical fibers Q , in superfluid he- 
lium Q, and in capillary waves Q, respectively. These 
discoveries indicate that rogue waves may be rather uni- 
versal. Certainly, we cannot harness the occurrence of 
rogue waves in the ocean due to their enormous destruc- 
tiveness. In optics, however, the optical rogue waves pro- 
duced in supercontinuum generation can be used to gen- 
erate highly energetic optical pulses @, fl0l - [T3} . 

Deep waves in the ocean and the wave propaga- 
tion in optical fibers can be described by the nonlinear 
Schrodinger (NLS) equation. Based on the model, the 
rogue wave phenomenon has been extensively studied, 
including rational solutions and their interactions 
l2l| . pulse splitting induced by higher-order modulation 
instability, and wave turbulence [22, [23| . A fundamen- 
tal analytical solution on the rogue waves is the Pere- 
grine solution (PS), which was first presented by Pere- 
grine 24]. PS is a localized solution in both time and 



space, and is a limiting case of Kuznetsov-Ma solitons 
[25l |26| and Akhmediev breathers [13] ■ Recently, the ex- 
citation conditions of PS have been demonstrated exper- 
imentally in optical fiber, and ex plic itly characterized its 
two-dimensional localization (28l. l29j| . It should be noted 
that the results are theoretically described by the NLS 
equation, which is valid for the picosecond pulses. When 
describing the characteristics of PS in the femtosecond 
regime, we must consider some higher-order effects, such 
as third-order dispersion (TOD), self-steepening and self- 
frequency shift, and so on. In this case, we should 
consider the higher-order nonlinear Schrodinger (HNLS) 
equation in the form (30j 
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where q(z,t) is the slowly varying envelope of the elec- 
tric field, z and t denote normalized propagation distance 
and the retarded temporal coordinate, respectively, and 
the parameters ai, 02, «3, 04, and 0:5 are the real con- 
stants related to the group velocity dispersion (GVD), 
the self- phase modulation (SPM), the third-order disper- 
sion (TOD), the self-steepening and the delayed nonlin- 
ear response effect, respectively. 

Generally, the HNLS equation ([T]) is not integrable. 
To solve Eq ([l}, we first consider the special parametric 
choice with «2 = 2/x 2 ai, 0:4 = 6/i 2 a3, 0:4 + 05 = so 



that Eq. ([T]) becomes [31 
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which is usually called the Hirota equation, where the 
parameter /i is a real constant. 

The equation was first presented by Hirota [3l| . and 
subsequently many researchers analyzed this equation 
from different points of view [33 - l39| . Recently, the rogue 
waves and the rational solutions of the Hirota equation 
have been discussed in the forms of the two lower-order 
solutions by employing the Darboux transformation tech- 
nique 40]. In this paper, based on the soliton solution on 
a continuous wave (CW) background for the integrable 
Hirota equation ©, we discuss the formation mecha- 
nism and the characteristics of the PS in the femtosecond 
regime. The results show that the PS exhibits a feature 
of temporal and spatial localization and can be excited 
by a small localized (single peak) perturbation pulse of 
the CW background, even for the nonintegrable HNLS 
equation ([T]). Based on this result, we discuss the dy- 
namics of the Peregrine rogue wave through a realistic 
example. 

The paper is organized as follows. In Sec. II we present 
the explicit process of the formation of the PS from the 
Kuznetsov-Ma solitons and Akhmediev breathers based 
on the soliton solution on a CW background for the Hi- 
rota equation and discuss the characteristics of the PS. 
Various initial excitations of the PS are discussed in Sec. 
III. Subsequently, in Sec. IV we investigate the influence 
of the self-frequency shift to the dynamics of the Pere- 
grine rogue wave by employing a realistic example. Our 
results are summarized in Sec. V. 



II. THE PEREGRINE SOLUTION INDUCED BY 
INTERACTION BETWEEN THE CONTINUOUS 
WAVE BACKGROUND AND SOLITON 

By employing a Darboux transformation, one can con- 
struct the soliton solution on a CW background for Eq. 

© as follows m m 



q(z,t) 



A + A 



+iA, 



a cosh 9 + cos tp 
cosh 9 + a cos uo 
b sinh 9 + c sin ip 



cosh4 



a cos if 



exp (iip c ) . (3) 



Here 



9 = MjT - ( Vl M R + u 2 Mi) Z, 
<p = M R T - {v 2 M R - i/xMj) Z, 
(f c =ujT+ [(2fi 2 A 2 - oj 2 ) ai + oj{6ij, 2 A 2 - uj 2 )a 3 ] Z, 

with T = t — to, Z = z — zq, and the coefficients are a — 
~2^ 2 AA S /^ 2 A 2 S + Ml), b = -2f,AM R /(^ 2 A 2 s + M%), 
and c = Mi / (ijlA s ) , which implies that Mi = as A s =0. 
And other coefficients v\ = [iA s [ai + (uj + 2ui s )a3], 
v-2 = (uj + u s )oti + (uj 2 + u)lj s + uj 2 - 2fj, 2 A 2 - fi 2 Al)a 3 , 
M R = {[((«. - uj) 2 + A^ 2 A 2 - n 2 A 2 s ) 2 + A^ 2 A 2 (oj s - 
uj) 2 ] 1 / 2 + ((uj s -uj) 2 +^ 2 A 2 - f i 2 A 2 s )} 1 / 2 /V2, and Mj = 
{[((«. - uj) 2 + A^ 2 A 2 - fi 2 A 2 ) 2 + ^ 2 A 2 s (uj s - uj) 2 ] 1 ' 2 - 
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FIG. 1: (Color online) The reduction process from the solu- 
tion © to the Peregrine solution Q as (A 3 ,ijj s ) — > (2A^,u>). 
(a") A s = A, uj a = 0.6a;; (b~) A 3 = 1A5A, uj s = 0.7a;; 
(c - ) A s = ISA, ljs = 0.8a;; (d _ ) A s = 1.99A, uj s = 0.85a;; 
(a+) A a = 3,4, cj s = 1.4a;; (b+) A s = 2.5/1, uj s = 1.3a;; (c+) 
As = 2.4/1, a; s = 1.2a;; (d+) A s = 2.01A, uj s = 1.15a;; and (e) 
the Peregrine solution given by ((4| with A s = 2 A and o; s = uj. 
The slope of red dashed lines is K v , and black solid lines is 



Kg. Here the parameters are A 
fi — 1, ai = 0.5, and 0:3 = 0.05. 



1, w = 1, to = 0, zq = 5, 



(K - uj) 2 + Ay 2 A 2 - ^A 2 )} 1 / 2 /^, and to, z , A, uj, 
A Sl and uj s are the arbitrary real constants, and with- 
out loss of generality we assume that A and A s are non- 
negative constants. The solution © includes two spe- 
cial cases. One is that as the amplitude A vanishes, it 
reduces to the solution q s (z,t) — A s e llPs sech#, where 
9 = [iA s T - nAsftujsCt! + (3uj 2 - fi 2 A 2 )a 3 ]Z and cp s = 
ip c + up = uj s T+ [(fi 2 A 2 s ~ uj 2 ) ai + uj s {3fi 2 A 2 s - uj 2 )a 3 ]Z, 
which describes a bright soliton solution with the max- 
imal amplitude A s . The other is that when the soliton 
amplitude A s vanishes, it reduces to the CW light so- 
lution q c (z,t) = Ae lVa . Therefore, in general, the exact 
solution q(z,t) in Eq. ([3]) describes a soliton solution 
embedded in a CW light background with the group ve- 
locity V sc = {viM R + v 2 Mi)/Mi [39]. Specifically, when 
a<3 = and taking ot\ = —1, /i = 1, the solution © co- 
incides with the result given in Ref . [4lJ , which is firstly 
derived from NLS equation by N. Akhmediev and V. I. 
Korneev [27|. 

In the limit (A s ,uj s ) ->• (2A,uj), Eq. ^ reduces to PS 
as follows 



q p (z,t) = Ae l 



4 + i8CZ 



1 + 4C 2 Z 2 + 4B 2 (T - DZ) 



(4) 



where B = [iA, C = 2B 2 [a x + 3uja 3 ) and D = 2uju\ 
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3 (a; 2 — 2£? 2 )a3, which is a rational fraction solution, and 
is first derived from the Kuznetsov-Ma breather (KMB) 
by Peregrine [23 - l26l | and so is called the Peregrine so- 
lution or the Peregrine rogue wave (Note that it is re- 
garded as the Peregrine soliton in Refs. [HI, Here 
(to, zo) presents the peak position of the PS. In order to 
understand the characteristics of the solution ©, Fig. 1 
presents the reduced process from the solutions ^ to the 
PS Q as (A s ,u) s ) (2A-,w) and (A s ,uj s ) ->■ (2A+,w), 
respectively. From them one can see directly that the 
solution ([3]) commonly exhibits breather characteristics, 
and the separation between adjacent peaks gradually in- 
creases as (A s ,uj s ) — >• (2A,lj), and eventually reduces 
into the PS. The most impressive feature on the PS is 
localized in both time and space. The peak position 
of the PS is fixed at spatiotemporal position z = zq 
and t — to during the reduction, and the peak power 
\q p (z ,t )\ 2 = (A + A s ) 2 9A 2 as A s 2A, which 
means that the PS with peak power 9 A 2 can be gener- 
ated by choosing an initial excitation properly. 

It should be pointed out that in the reduction process 
mentioned above, the physical mechanism for the forma- 
tion of the PS has some differences. Figures l(a~)-l(e) 
demonstrate a localized process of CW background along 
the slope direction K v = Mij/(j/ 2 M^— v\ Mi), while Figs. 
l(a + )-l(e) show a periodization process of bright soliton 
along the slope direction Kg = Mj / (i/iMr + ^Mi). Fur- 
thermore, Fig. 2 presents the contour plots of Kg and K v 
as a function of A s and lo s for given A and lj. Note that 
K v is infinite as v^Mji — v\Mj — except for the point 
(2A,u>), which corresponds to the line t = to and appears 
at the sixth curve in Fig. 2(b). When (A s ,uj s ) — »• (2A,ui), 
the limitations of Km and Kg do not exist, and at this 
point the solution is localized along both slope directions 
K v and Kg, namely, the PS appears. Therefore, the PS 
should be a middle state in the process of the localiza- 
tion of CW converting into the periodization of the bright 
soliton, as shown in Ref. [17j . 

In order to understand better the reduction processes 
of the PS, we consider a special case for the solution 
(|3|) with uj s — lj. In this case, the solution (j3|) has two 
different expressions. When A 2 < AA 2 , the solution ^ 
can be written as 



q(z,t) = fin 



f2i{ cosh 9 + iA s sinh 9 
2 A cosh 9 — A s cos op 



-A 



(5) 



where 9 = /i 2 A s il R (ai + 3u>a 3 )Z and ip — [&IrT — 
/j,Q R [2ujai + (3w 2 - 2\i 2 A 2 - fi 2 A 2 )a 3 ]Z with £l R = 
a/ 4A 2 — A 2 S being a modulation frequency. Here Kg = 
wAK v = l/[2oj ai + (3uj 2 -2fi 2 A 2 - n 2 A 2 s )a 3 ]. Therefore 
the solution ([5]) is periodic with period t pcl — 2tt/([iHr) 
along the t axis and localized along the z axis, and is 
usually called an Akhmediev breather (AB), which can 
be considered as a modulation instability process [42[ . In 
this case, the solution ([5]) can reduce to the CW back- 
ground q c {z, t) as A s — > 0, and with the increasing of A s 
the CW background is gradually localized due to the in- 
teraction with the soliton and forms a periodic breather 





0.50 


0.20 






' o >z — 


" ).02O T ^S^M 




— 'ojp. -"' 7 




0.40 1 




0.30 / / 


0.60 

(a) N 



14 > 30 I 

?i.e\ 

\ \ \ I 


-0 90 

-2.5 

/ / -0.40 
\// -° A0 " 


3.o' / 


I \\ -0.20 

-0.90 
15 \ \> 



2.0 
A 



2 

A 



FIG. 2: (Color online) The contour plots of (a) Kg and (b) K v 
as a function of A 3 and u; s , where the parameters are A = 1, 
w = 1, jLt = 1, ai = 0.5, and ay = 0.05, and the contour 
values (from left to right) are 0.02, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 
0.7, 0.9, 1.1, 1.3 in (a), and the contour values are 1.1, 1.2, 
1.4, 1.8, 3, -2.5, -0.9, -0.4, -0.2, -0.1 in (b), respectively. 



with period i por , i- e, the AB [see Fig. 3(a)], eventually 
when A s tends to 2A, the solution ([5]) becomes the PS. 
This process represents the reduction of CW— >AB^PS. 
Also, it can be described by the ratio between the period 
t pel and the temporal width 5t as follows: 



t 



per 

Sto~ 



2tt 



cos 



As _ 9 ^4_ 

A Z A S 



(6) 



where the temporal width Sto is defined as the width 
from zero-value of intensity to adjacent peak 28] . Fig- 
ure 3(b) presents the dependence of the ratio t per /5to on 
A s j A. From it, one can see that when A s /A approaches 
2, the ratio t pey ;/Sto tends to infinity, which implies that 
the separation between peaks is more and more large, 
eventually resulting in the localization along the t axis, 
and forms the PS. From the expression ©, one can see 
that this process does not depend on the higher-order 
parameter a 3 , which means that when the higher-order 
effects with the conditions «2 = 2/i 2 ai, = 6fi 2 a 3 , and 
Q4 + as =0 simultaneously appear in the optical fiber, 
they do not influence the characteristics of the PS. It 
should be emphasized that the Peregrine solution gener- 
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FIG. 3: (Color online) (a) and (c) The evolution plots of 
the solutions ([5]) and (0, respectively, where the parameters 
are the same as in Fig. 1 except for A s = 1.5A in (a) and 
A s — 2.5A in (c); (b) and (d) the dependences of the ratio 
tpcr/Sto and z peT /Szo on A s /A, respectively. 



ated by the process of the reduction of CW— !>AB— >PS has 
been already studied theoretically and experimentally in 
the framework of the NLS equation [lj], l26U29j . Here 
we presented the corresponding descriptions in the fem- 
tosecond regime. 

When A 2 S > AA 2 , the solution ([3]) can reduce to the 
following form 

y ' \ A s cosh 9- 2A cos tp J w 

where 9 = ^VtiT '- ^/[2o;ai + (3u; 2 -2[i 2 A 2 - [i 2 A 2 )a 3 ]Z 
and <p = fi 2 A s n I (a 1 + 3uja 3 )Z with 0/ = A 2 - AA 2 . 
Here Kg = l/[2ua x + (3w 2 - 2fi 2 A 2 - ^ 2 A 2 s )a 3 ] and 
K v = 0. Especially since A s ^> 2 A, namely f2/ ~ A s , the 
solution (J7J can be approximated as q{z,t) w q s (z,t) — 
q c (z, t), which is the superposition of a CW solution and 
a bright soliton with the larger amplitude A s . From the 
solution (O, one can see that it is a periodic function 
with period z pcl — 2ir/ [i 2 A s £li(ai + 30-10:3) along the z 
axis and is localized along the t axis, possessing the pe- 
riodic peaking property of the field amp litude like the 
Kuznetsov-Ma soliton (KMS) [H [H, Elf , and so is usu- 
ally called Kuznetsov-Ma soliton, as shown in Fig. 3(c). 
In this case, the solution (|7|) can reduce to the bright 
soliton q s (z,i) as A — > 0, and with the increasing of A, 
the bright soliton is periodized due to the interaction 
with the CW background and forms a KMS, and eventu- 
ally reduces into the PS. This process represents the re- 
duction of bright soliton— ^KMS—s-PS. Similarly, this pro- 
cess can be described by the ratio between the period 
z pcr = 2ir J [i 2 A s Q,j{ai + 3^0:3) and 5zq in the form 



^pcr 

Sz 



cos 



A S (-7A 2 -2AA B +A 2 S ) 
2A(A 2 -2AA S -A 2 S ) 



(8) 



where 5zo is defined as the distance between two cor- 
responding locations of half of the peak intensity along 
the slope direction Kg. Figure 3(d) presents the depen- 
dence of the ratio z peT /5zo on A s /A. From it, one can see 
that when A s /A approaches 2, the ratio z poT /Szo tends 
to infinity, which means that the distance between peaks 
is more and more large, eventually resulting in the for- 
mation of the PS. Similarly, the expression JSJ does not 
depend on the parameter a 3 . 

Another feature of PS can be described by the z- 
independent integral and the energy exchange between 
the PS and the CW background. Indeed, by introducing 
the light intensity against the CW background as follows, 



E c (t, z) = \q p (t, z)\ 2 - \q p (±00, z)\ 2 



A 



32C 2 Z 2 - 32B 2 (T - DZy 
[1 + 4C* 2 Z 2 + 4B 2 (T - DZ) 2 ) 2 ' 



(9) 



it can be shown that it possesses the ^-independent in- 
tegral property, i. e. f_ °° E c (t, z)dt — 0. From the con- 
dition E c (to ± \/2B,zq) — 0, one can define the width 
of PS as l/B, and have integrals f^j^g E c (t, z )dt = 

4A 2 /B and f:^ 2B E c (t, z Q )dt + f™ /2B E c (t, z Q )dt = 

— AA 2 /B. These results show that the energy of PS with 
stronger intensity mainly concentrates in its central in- 
terval (— 1/2B, 1/2B), but as a result of z-independent 
integral properties, it loses the same energy in the back- 
ground so that the relation 51 + 5*2 = S3 holds, as shown 
in Figs. 4(a)-4(c). It should be noted that because 
B = fiA is independent of the higher-order terms in Eq. 
@, the higher-order effects do not influence this property 
on the PS. 

The energy exchange between the PS and the CW 
background is of the form 



E e (z) 



+ OO 



\q p (t,z) - q p (±oo,2)| dt 
4A 2 tt 

By/1 + 4[2S 2 (ai + 3oja 3 )} 2 Z 2 ' 



(10) 



From the expression (fT0|) . one can see that E e {z) is a 
periodic function of z, which differs from the periodic 
exchange between the bright soliton and the CW back- 
ground in Eq. ([3]) 39]. Also, one finds that E e (z) is 
monotonously increasing as z < zq and is monotonously 
decreasing as z > zq; thus it takes a maximal value at 
z — zo, as shown in Fig. 4(d). This result shows that 
at z — zq the energy exchange between the PS and the 
CW background reaches maximum. In order to under- 
stand the influence of the higher-order effects on the PS, 
Fig. 4(d) presents the evolution plots of E e (z) for differ- 
ent TOD parameter a 3 . From it, one can see that the 
increase of the TOD parameters a 3 can enhance the rate 
of the energy exchange, as shown in Figs. 4(a) and 4(b). 



5 





Time/ Distance z Timet Time/ 



FIG. 4: (Color online) The evolution plots of the Peregrine 
solution given by Eq. @ for 03 = in (a) and Q3 = 0.2 in 
(b), respectively, (c) The distribution of the light intensity 
against the CW background given by Eq. ((9J at z — zq. (d) 
The evolutions of the energy exchange between the PS and 
the CW background for the different TOD parameters. Here 
the parameters are A — 1, ui = 1, to = 0, zo = 5, p = 1, and 
on = 0.5. 



III. INITIAL EXCITATIONS OF THE 
PEREGRINE ROGUE WAVES 

In this section we will discuss the initial excitations 
of the Peregrine rogue wave. We start with considering 
the excitation of the Peregrine rogue wave based on the 
solutions ([5]) and ([7]). For solution (|SJ), by linearizing its 
initial expression, one finds that the initial expression of 
the solution (|5|) can be approximated by 

q±(0,t)tt(p± + e ± CO S cp)e i ^ > t \ (11) 

where p± = (2 A 2 - A 2 S =p iA s VL R )/(2A) with \p±\ = 
A, e± = A s n R {Q R T iA s )/(2A 2 )exp[Tp 2 A s (a 1 + 
3wa3)f2_RZo]i an d <p = p^rT + pQ R [2uja>i + (3w 2 — 
2p 2 A 2 — p 2 A 2 s )az\z . Here the subscripts "+" and "— " 
correspond to the case of a\ + Suias > and ol\ + iuia^ < 
0, respectively. It can be shown that e± — > as A s — >• 2A. 
Thus the expression (jTTJ) can be regarded as an initial 
condition with a small periodic perturbation of back- 
ground with the period i por . The numerical simulations 
show that the evolution of the exact solution ([5]) can be 
well described by the initial approximation (ITTj) (also see 
Ref. [39j]). Here one makes use of the initial approxima- 
tion (jlip to excite the Peregrine rogue wave as A s closes 
to 2A. 

Figure 5 presents the evolution plots of the numerical 
solution of Eq. with the initial condition (jTTJ) and 
the comparisons of the intensity profile between numeri- 
cal and exact Peregrine rogue waves at peak position for 
a 3 = and 113 ^ 0, respectively. From it, one can see 
that when A s approaches to 2A, the initial condition pT|) 
evolves into a string of near-ideal Peregrine rogue waves. 



FIG. 5: (Color online) The evolution plots of the numerical 
solution of Eq. ((2| with the initial condition and the 
comparisons of the intensity profiles of numerical and exact 
results at peak position, respectively: (a,c) 03 = 0, (b,d) 
«3 = 0.2. Note that in (c) and (d), the blue dotted curves 
are the exact results given by Eq. at z = Zn and the 
black solid curves are the numerical results at z — 5.9 in (a) 
and z = 5.2 in (b) [the white dotted lines in (a) and (b)], 
respectively. Here A„ = 1.95, and other parameters are the 
same as in Fig. 4. 



Although the theoretical results reveal that there is a 
Peregrine rogue wave only in the limitation of A s —> 2A, 
in practice, this cannot be implemented due to e± = 
as A s = 2 A. Furthermore, from Figs. 5(a) and 5(b), it 
can be seen that the evolution for long distances shows 
the breakup of the Peregrine rogue wave, which implies 
that the Peregrine rogue wave is unstable. Also, from 
Figs. 5(c) and 5(d), one can see that the presence of the 
higher-order effects did not markedly influence the in- 
tensity distribution of the Peregrine rogue wave at peak 
position, except for a displacement of the peak position, 
but shortens the distance of energy exchange, as shown 
in Figs. 5(a) and 5(b). This is in agreement with that 
suggested in Fig. 4(d). 

In the following, we discuss the initial excitation in- 
duced by the solution 0. In this case, one find that the 
solution can take the following particular form 

q(zi,t) = (-A + iQ/sech^e^* 1 '*) 

= Ae*^* 1 +ft / sech(9e l M z i< t )+f] (12) 

at the location z\ = (2nir + ir/2)/[p 2 A s rii(ai +3wa3)] + 
z , n = 0,±1,±2, - • •, where = /ifi/T- [2wai + (3w 2 - 
2p 2 A 2 - p 2 A 2 s )a 3 ](2mT + n/2)/[pA s (a 1 + 3wa 3 )]. With- 
out loss of generality, here we take n = 0. The expres- 
sion (fl~2"j) is the superposition of a CW and an aperiodic 
hyperbolic secant function, especially as A s — > 2A, ilj 
tends to zero. This means that the expression ([T2j) can 
be regarded as an initial condition with a small aperiodic 
(simple peak) perturbation of background, which differs 
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FIG. 6: (Color online) The evolution plots of the numerical 
solution of Eq. ((2|) with the initial condition (|12[) and the 
comparisons of the corresponding intensity profile at peak 
position, respectively: (a,c) a-j — 0, (b) and (d) Q3 = 0.2. 
Note that in (c,d), the blue dotted curves are the exact re- 
sults given by Eq. @ at z — zo and the black solid curves 
are the numerical results at z = 10.2 in (a) and z = 6.4 in 
(b) [the white dotted lines in (a) and (b)], respectively. Here 
A s = 2.05, to = 0.5, and other parameters are the same as in 
Fig. 4. 
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FIG. 7: (Color online) The numerical evolution plots of the 
initial condition 1)131) : (a) for NLS equation with 03 = 0; (b) 
for Hirota equation ((2)1 with «3 = 0.1; (c) for HNLS equation 
Q with ay = 0.05, a 4 = 0.285, and a 5 = -0.315; and (d) 
the distributions of the intensity profile at z — 6.975 in (a), 
z = 6.825 in (b), and 2 = 6.9 in (c). Here ai = 0.5, A = 1, 
e = 0.05, and cr = 0.05. 



type perturbation pulse as follows 



from the superposition of a CW and a periodic pertur- 
bation in the expression Here we make use of the 
expression (|12l) as an initial condition to investigate the 
excitation of the Peregrine rogue wave. Similarly, the nu- 
merical simulations show that the evolution of the exact 
solution © can be well described by the initial condition 
(fT2j) except for a translation. 

Figure 6 shows the evolutions of the numerical solution 
of Eq. with the initial condition (fTIZj) and the com- 
parisons of the intensity profile of numerical and exact 
Peregrine rogue waves at peak position for a?, = and 
a 3 7^ 0, respectively. From it, one can see that the initial 
condition p^|) evolves into a near-ideal Peregrine rogue 
wave, which differs from that shown in Figs. 5(a) and 
5(b). Similarly, Figs. 6(a) and 6(b) show the breakup of 
the Peregrine rogue wave for long distance, which implies 
that the Peregrine rogue wave is unstable. 

Comparing results in Figs. 5 and 6, we find that a small 
localized (single-peak) perturbation pulse of CW back- 
ground can excite a Peregrine rogue wave. Indeed, every 
peak in the initial expression can excite a Peregrine 
rogue wave; thus a periodic perturbation of background 
resulted in the generation of a string of the Peregrine 
rogue waves, as shown in Figs. 5(a) and 5(b). So we can 
suggest that the Peregrine rogue wave can be excited 
by the interaction between the CW and a small local- 
ized (single-peak) perturbation pulse. As an example, 
we consider a simple initial condition with a Gaussian- 



q(0,t) = A + ecxp(-crt 2 ), (13) 

where e is a modulation amplitude and a small quantity. 
The numerical simulations show that when the width of 
the initial perturbation pulse is wide enough, i. e, the pa- 
rameter a is small enough, the Peregrine rogue wave can 
be excited by the initial condition (1131) even for a noninte- 
grable HNLS equation fTJ). Figure 7 shows the evolution 
plots of the numerical solution for Eq. (JTJ) with the initial 
condition (|13p . which includes a numerical evolution of 
the initial condition (| 13[) under the nonintegrable case, 
as shown in Fig. 7(c). From Fig. 7(d), one can see that 
the main characteristics of the Peregrine rogue wave are 
maintained. Certainly, the Peregrine rogue wave is un- 
stable. This result can be used to understand extreme 
localized events in ocean. 



IV. THE INFLUENCE OF THE 
SELF-FREQUENCY SHIFT TO THE DYNAMICS 
OF THE PEREGRINE ROGUE WAVE 

It should be noted that the equation (JTJ) does not in- 
clude the self-frequency shift effect arising from stimu- 
lated Raman scattering because the parameter a§ is a 
real number. In this section we discuss the influence of 
the self-frequency shift effect on the dynamics of the Pere- 
grine rogue wave. In this case, the model governing pulse 
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propagation can be written as [4c 
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where r and £ represent the temporal coordinate and 
the propagation distance, /?2 is the group velocity disper- 
sion, /?3 is the third-order dispersion, 7 is the nonlinear 
coefficient of the fiber, and Tr is the Raman time con- 
stant. By introducing the dimensionless transformations 
Q(£,r) = y/P^q(z,t), t = t/T , and z = £/L D with the 
dispersion length Lo = Tg/|/?2|, Eq. (| 14[) becomes the 
form of Eq. ([T]) with the coefficients a.\ — — ^2/(2 |/?a|), 
a 2 = lL D Po, «3 = Aj£d/(6T 3 ), a 4 = -7L D P /(w T ), 
and «5 = —{jTrLdPo/Tq, respectively. It should be 
pointed out that here as is a complex number which de- 
scribes the self-frequency shift effect arising from stimu- 
lated Raman scattering. 

As an example, here we use realistic parameters for a 
highly nonlinear fiber at 1550nm with the group velocity 
dispersion fa = —8.85 x 1CP 1 ps 2 /km, the third-order 
dispersion fa = 1.331 x 1CP 2 ps 3 /km, and the nonlinear 
parameter 7 = 10 W~ 1 -km _1 [28|]. Thus, for a given 
initial power P , the parameters ai, 0:2, CI3, 014 and a$ 
can be determined. Note that in our simulations we take 
a 2 = 1 by choosing a temporal scale T = [\fa \/ (jPq)] 1 ^ 2 , 
and the Raman time constant Tr = 5fs when the self- 
frequency shift is considered. We still take Eq. ([T3"l) as 
the initial condition, in which the realistic width of the 
initial perturbation pulse is T\ = To/^/2a. 

Figure 8 presents the numerical evolution plots of the 
initial condition (|13[) with the initial perturbation pulse 
width Ti = 2ps for the different initial power P . From 
Fig. 8 it can be seen that, in the absence of the self- 
frequency shift effect (Tr = 0), the Peregrine rogue wave 
in turn splits into two subpulses, three subpulses, and so 
on, and can split into more subpulses for higher initial 
power, as shown in Figs. 8(a) and 8(b). These results 
are similar to the pulse splitting induced by higher-order 
modulation instability based on NLS equation in Ref. 
22], but here the complex splitting dynamic evolutions 
are excited by a linear superposition of the CW and a 
small localized (single peak) perturbation pulse, and the 
higher-order effects, such as the third-order dispersion 
and the self-steepening, are included. However, when 
the self-frequency shift effect is considered, such splitting 
of the Peregrine rogue wave no longer appears, as shown 
in Figs. 8(c) and 8(d). It is surprising that, in this case, 
the dynamics of the Peregrine rogue wave mainly exhibit 
a peak changing propagation characteristic of the field 
amplitude except for some of the small radiations, and 
has an increase in peaking value, as shown in Figs. 8(e) 
and 8(f). This property can be used for generation of the 
higher peak power pulse. 

Furthermore, the dependences of the peak position 
£pcak of the excited Peregrine rogue wave on the mod- 
ulation amplitude e are considered, as shown in Fig. 9. 
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FIG. 8: (Color online) The numerical evolution plots of the 
initial condition (|13[l for the different initial power: (a,b) the 
absence of self-frequency shift (Tr = 0), (c,d) the presence of 
self-frequency shift (Tr = 5 fs), and (e,f) the corresponding 
power density distributions at the different distance labeled 
by the white dotted lines in (c) and (d), respectively. In (a), 
(c), and (e), the initial power Po — 0.3W, and in (b), (d), and 
(f), the initial power P = 0.4W. Here e = 0.2 and A = 1. 




FIG. 9: (Color online) The dependences of the peak position 
Cpcak of the excited Peregrine rogue wave on the modulation 
amplitude e for (a) the different initial power and a given ini- 
tial perturbation pulse width T\ — 2ps, and (b) the different 
initial perturbation pulse width T\ and a given initial power 
Po — 0.4W. Here the solid, dashed and dotted curves corre- 
spond to the results of Tr = 0, and the open diamonds, circles 
and triangles correpond to the results of Tr = 5fs, and the 
parameter A = 1. 



From it one can see that for a given initial perturba- 
tion pulse width T\ or initial power Pq, the peak position 
Cpcak of the excited Peregrine rogue wave is a decreasing 
function of the modulation amplitude e. Thus one can 
control the position of the excited Peregrine rogue wave 
by suitably choosing the initial perturbation pulse width 
or the initial power. Also, we find that the self-frequency 
shift effect does not influence the peak position £ pca k of 
the excited Peregrine rogue wave. 
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V. CONCLUSIONS 

In summary, based on the soliton solution on a CW 
background for an integrable Hirota equation, we have 
presented the reduction mechanism and the main char- 
acteristics of the Peregrine rogue wave in the propaga- 
tion of femtosecond pulses of optical fiber. The results 
have shown that there exist two processes of the forma- 
tion of the Peregrine rogue wave: one is the reduction 
of CW^AB^PS, which is the localized process of the 
CW background, and the other is the reduction of bright 
soliton— >KMS->PS, which is the reduction process of the 
periodization of the bright soliton. The characteristics of 
the Peregrine rogue wave have been exhibited by strong 
temporal and spatial localization. Also, the initial exci- 
tations of the Peregrine rogue wave have been discussed. 
The results have shown that a Peregrine rogue wave can 
be excited by a small localized (single peak) perturba- 



tion pulse of CW background, even for the nonintegrable 
HNLS equation. This means that the Peregrine rogue 
wave is a result of interaction between the continuous 
wave background and soliton. Furthermore, the numeri- 
cal simulations have shown that the Peregrine rogue wave 
is unstable. Finally, through the study for a realistic 
highly nonlinear fiber, it has been found that the self- 
frequency shift influences the dynamics of the Peregrine 
rogue wave. The results show that in the absence of the 
self-frequency shift, the Peregrine rogue wave can splits 
into several subpulses; however, when the self-frequency 
shift is considered, the Peregrine rogue wave no longer 
splits and exhibits mainly a peak changing and increas- 
ing propagation property. 

This research was supported by the National Natural 
Science Foundation of China under Grant No. 61078079 
and the Shanxi Scholarship Council of China under Grant 
No.2011-010. 
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